The origin of dynamic scaling in dilute polymer solutions 
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The hydrodynamic radius of a polymer chain, obtained using Brownian dynamics simulations of the contin- 
uum Edwards model, is found to obey a crossover in the excluded volume parameter z, which is significantly 
different from that observed for the radius of gyration. It is shown that this difference arises from contributions 
due to dynamic correlations to the diffusivity, which are ignored in the commonly used definition of hydrody- 
namic radius based on the Kirkwood expression. The swelling of the hydrodynamic radius from the (9-state, 
obtained from simulations, shows remarkable agreement with experimental measurements. 
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The mean size of a polymer in a dilute solution can be ob- 
tained by static measurements that yield the radius of gyration 
i? g or by dynamic experiments that provide the diffusion co- 
efficient D, An outstanding problem yet unresolved is the 
observed difference in the scaling with molecular weight be- 
tween static and dynamic measurements UJ. For instance, in 
the limit of large molecular weight M, R g scales as M 59 , 
while the hydrodynamic radius i?ii(oc D~ r ) scales as M 57 
Q|. In addition, in the crossover regime between the 0-state 
and the good solvent limit, the growth of R g and i?u rel- 
ative to the 0-state, in terms of the scaling variable z — 
vq (1 — Tg/T) \fM, is observed to be significantly different. 
Here, vq is a chemistry dependent prefactor, T is the temper- 
ature, and Tg is the temperature at the 0-state. Theoretical 
attempts to explain the origin of this anamalous dynamic scal- 
ing behavior have mainly based the definition of i?n on the 
Kirkwood expression for diffusivity |3J |4|, |5J, |a] . In this work, 
we show using exact Brownian Dynamics simulations of the 
continuum Edwards model [7], that (i) the use of Kirkwood's 
expression in fact leads to nearly identical scaling in the two 
cases, (ii) in order to exhibit dynamic scaling that is distinct 
from static scaling, i?H should be based on the long-time dif- 
fusivity, and (iii) the latter definition of leads to excellent 
agreement with experimental observations. 

Theoretical treatments of the static and dynamic properties 
of polymers in a dilute solution are typically carried out with 
a bead-spring chain model, which consists of N beads con- 
nected by linear springs. The swelling of a polymer in a good 
solvent is usually represented in terms of the ratio of the size 
in the presence of excluded volume (EV) interactions between 
the beads to size in the #-state. The swelling of the radius of 
gyration is denoted as, a g = R g /R^ and that of the hydrody- 
namic radius as an = Rh/R^. The hydrodynamic radius 
is often approximated by a static measure constructed from the 
Kirkwood estimate of the diffusivity |8] 
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Here, 



is the distance between the beads /i and v. It is 



where, R\ is the inverse radius (or sometimes known as the 



observed from experiments that an < a g , and for the typ- 
ical maximum molecular weights studied, the effective scal- 
ing exponent (defined as v cS = 1/2 + 5 log a/d log M) is 
jy cff = 0.59 for a g and v cS — 0.57 for an- An early expla- 
nation of this anomalous dynamic scaling behavior employed 
a model of self-excluding Gaussian blobs 0] to show that the 
segments within the Gaussian blob have more influence on 
l/Ri than on R g . It was shown that this leads to a swelling in 
Ri, defined as a.\ — R\jR\, which obeys a.\ < a g . The effec- 
tive exponent for Ri was shown to reach the asymptotic value 
of v = 0.59 much more slowly than that for R g , and an ap- 
proximate estimate provided v cS = 0.57 for Ri for the exper- 
imental system considered. This idea and other modified blob 
theories have subsequently been discounted by Monte Carlo 
simulations of excluded volume chains which show 

that the model ignores the excess swelling of internal parts of 
the chain and that in fact the effective exponent is significantly 
closer to 0.59 than 0.57 full , a fact supported by renormali- 
sation group (RG) calculations 1 3] and our own results shown 
here. 

The crossover scaling behaviour (from the 6>-limit to the 
good-solvent limit) of static properties, such as a g is a well 
studied problem. One of the most successful and widely 
used approaches is the Edwards continuous chain (N — > oo) 
representation of the bead-spring model iTJ]. RG calcula- 
tions have shown that this model correctly captures the static 
scaling behavior for temperatures asymptotically near the 9- 
temperature, in close agreement with experimental observa- 
tions HGilllUll]. However, RG calculations Q of the 
crossover of dynamic swelling ratios based on the definition 
of the Kirkwood diffusivity lead to a.\ w a g (to within 1%). 

An alternative explanation for the lower value of cx\ is the 
existence of a draining effect in hydrodynamic interactions 
01 . This implies the introduction of another (draining) pa- 
rameter into the theory. However, a recent compilation of dy- 
j i 

namic data 1 15] shows a significant degree of universality in 
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the crossover in z, clearly indicating the absence of another 
parameter (or that the chains are non-free draining). 

In this work, we have reexamined the dynamic crossover 
problem, (restricting our attention to the Edwards' continuum 
chain description) in the light of three recent developments. 
First is the introduction of a scheme to simulate the continuum 
Edwards model using the Brownian dynamics (BD) method 
lllfill . which permits the study of dynamical and rheological 
properties. In this scheme, a repulsive Gaussian EV poten- 
tial is used to mimic the <5-function potential in the Edwards 
model, and the mapping of the model parameters to the EV 
parameter z is shown to occur in the continuous chain limit 
11711 . The second is the advance made in the treatment of 
fluctuating HI using fast approximation methods 11811 . The 
third is the resolution of a long standing problem with regard 
to the diffusivity of a polymer chain. Liu and Diinweg Il9tl 
have numerically proven that for chains with fluctuating HI, 
the long time diffusivity is lower than the short time Kirk- 
wood estimate Dk, and is given by D = Dk — D\, where D\ 
is the equilibrium time correlation function of the drift term 
in the stochastic differential equation governing the motion of 
the beads of the chain [20]. These developments permit us to 
compute the diffusivity from the Edwards model at a finite z 
and obtain the crossover exactly in the non-free draining limit. 
We show that in the crossover regime a.\ w a g , and an ob- 
tained from the long-time diffusivity is able to reproduce the 
observed anomaly in the swelling of hydrodynamic radius. 

To simulate the continuum Edwards model, we begin with a 
discrete bead-spring chain with N beads subject to Brownian 
motion. The dimensionless displacement of a bead position 
AR^i during a finite time increment At is taken to be |21]: 

AR^ = - F v j At+-^= B Muij AW vj , (3) 

where, the summation convention is implied for repeated in- 
dices, and a Greek subscript [i,v.... = [1,2, ... ,N] denotes 
the bead index and a Roman subscript = [1,2,3] de- 

notes a Cartesian coordinate. Here, is the diffusion 

tensor, F^, is a matrix of the total body force acting on the 
beads, B^ij is taken as the square root matrix of de- 
fined as B^QikBvQjk = D^vij, and AW a % is an increment to 
the Weiner process W V j. In the absence of EV interactions, 
the distribution of bead connector vectors is assumed to be 
Gaussian. A repulsive narrow-Gaussian potential is used to 
regularise the <5-function EV potential 1 17] considered in the 
Edwards model |7]. The potential is parameterised by z* rep- 
resenting the strength and d* representing its width llal . The 
solvent quality (EV parameter) is given by z = z* \ ~N, ob- 
tained in the simultaneous limit of N — > oo and z* — > 0. 
In Ref. llal the expression d* = K z* 1 ^ 5 (where if is a con- 
stant) was found to be a computationally efficient way to reach 
the limit d* — ► 0, which recovers the (5-function. This method 
was shown to accurately predict the swelling of the gyration 
radius 1 16] observed in experiments [22]. To simulate hydro- 
dynamic interactions (HI) between the beads we have used the 
usual Rotne-Prager-Yamakawa expression for the HI tensor, 



and an accelarated method of finding the square root matrix 
Bavin dl. 

These and other details of the simulation are dis- 
cussed elsewhere 12111 . 

Fixman i2oll established from a preaveraged treatment for 
Gaussian chains that dynamic correlations play a minor role 
in the diffusivity of chains. He showed that the long time dif- 
fusivity is given by D = Dk — D±, where the contribution 
Di > is related to the time-correlation function of the drift 
term in Equation 0, and is only about 2% of Dk- Defining 
A i = (X^ D »»ij Fvj)/^, we have 

OO 

Dl = 3^ J dt (MVMt))- (4) 

o 

Liu and Diinweg 1 19], who recently performed simulations to 
obtain diffusivity D directly, and through the Fixman's rela- 
tionship, effectively proved the relation by showing numeri- 
cally (for excluded volume chains) the variation of diffusivity 
from its short time value Dk to the long time limit. Even 
in this case, it was found that the contribution of D\ is only 
a small fraction of D. This is one of the prime reasons for 
the neglect of dynamic correlations in the diffusivity and the 
use instead, of the Kirkwood estimate Dk — an approxima- 
tion which we show below has gone unnoticed in attempts to 
explain dynamic scaling. 

We have validated Fixman's expression in 0-conditions 
with fluctuating HI, by carrying out similar simulations to that 
in Ref. 1 19]. In this work, we have used the computationally 
easier method of evaluating D, which is to evaluate the static 
average Dk and obtain D\ from the time-correlation function. 
This immediately enables one to study the crossover of both 
the static estimate and the dynamically correlated value. 

We define two swelling ratios for the two different size 
measures obtained from the diffusivity. The commonly used 
Kirkwood estimate gives the swelling of the inverse radius 
ol\ = Ri/Rf = D k /Dk, in which the second equation is 
obtained in the limit N — > oo (which is true for the con- 
tinuum model). The hydrodynamic radius obtained from the 
long-time diffusivity has a swelling ratio qh = D e /D. ot\ 
is obtained accurately from free-draining runs by carrying out 
simulations at two different values of K = 1.0 and 1.5. qh is 
obtained from separate runs incorporating fluctuating HI and 
is carried out for two values of HI parameter (which is the 
dimensionless radius of a bead), h* — 0.19 and 0.30. The re- 
sults from simulations performed at finite TV are extrapolated 
to N — > oo, keeping z and h* constant. This limit corresponds 
to the continuum Edwards model in the non-free draining limit 
(h* VN -> oo). 

The universal nature of the swelling ratios ctj and an is 
seen in Figure^ The independence of the extrapolated value 
of ai from the width of the Gaussian repulsive potential indi- 
cates z is the only relavant parameter. A similar behaviour is 
observed for an (not shown here). To calculate an we find it 
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FIG. 1: Universality at a fixed solvent quality z in the continuous 
chain limit N — > oo. (a) The independence of the extrapolated value 
of ai from the width of the Gaussian repulsive potential and (b) A 
similar independence of (7rd from the HI parameter ft* at a solvent 
quality z = 1. 
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FIG. 2: Universal solvent quality crossover of swelling of two mea- 
sures of polymer size based on diffusivity (inverse radius ai and hy- 
drodynamic radius an) compared with the crossover of swelling of 
the radius of gyration (a g ). 



convenient to use the following expression 
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where, {/rd = Rg/Rn- This separates out the weak depen- 
dence of an on h* at finite TV into the two strong dependences 
of each of J7rd and U^ D . a g is independent of h* and can be 
obtained accurately from free-draining simulations along with 
a\. Our simulations of Gaussian chains with fluctuating HI 
yield a value of [7p D = 1.38 ± 0.01. a g is determined from 
free-draining simulations Jlrill . C/rd is determined by carry- 
ing out simulations at constant h*, and the extrapolations to 
N — > oo are carried out for two values of h* = 0.19 and 0.30 
as shown in Figure^). The two extrapolations can be seen to 
approach a common limiting value, which are not as close to 
each other as observed for the EV parameter. We believe that 
this is a limitation of the maximum N that we have consid- 
ered in the extrapolation, and a higher N (which is difficult 
to simulate presently) should lower the margin of difference 
in the extrapolated value. In the results to follow, we indicate 
this margin through an error-bar in the data points. 

The complete crossover of the two swelling ratios is shown 
in Figure |3 and is compared with that of a g . It is observed 
that aj w a g within error bounds for a large range of z, i.e., 
the cross-over of static size measures is nearly identical. This 
is a new result, calculated exactly for the first time for the 
continuum Edwards model. This is in contrast to the simple 
scaling arguments of the thermal blob model |4]. It is clear 
that the "spatial-crossover" explanation of the blob model is 
insufficient to show the observed difference in the crossover 
of the swelling of hydrodynamic radius. On the other hand, 
we see from Figure [2] that the crossover obtained from the 
long-time diffusivity is significantly different and o?h < ol s 
for all values of z. Recognition of this fact is important since 
most of the theoretical calculations are based on the static size 
measure (-Dk) ignoring dynamic correlations (D{). 

The theoretical crossover behaviour can be compared with 
experiments in two ways. In one method we can compare 
the z-crossover of the swelling ratios, for which it is required 
to estimate the phenomenological parameter z for the exper- 
imental system, which has a chemistry dependent prefactor. 
In the alternative approach, we can directly compare with two 
experimentally measured quantities, without the need to esti- 
mate z. We have chosen the most recent data (to our knowl- 
edge) on the hydrodynamic radius, collated in Ref. 1131 . in 
terms of a solvent quality they denote as z. We found that by 
taking z = 0.85 z, we can shift the experimental data for the 
swelling of radius of gyration a g on to our theoretical curve. 
We then compare the crossover of an using the same shift 
factor. Figure [5] shows this comparison leads to remarkable 
agreement. It is also observed that the crossover of aj lies 
outside the region of experimental accuracy. This provides a 
substantial evidence that dynamic correlations are important 
to describe the observed crossover behavior. We have also 
made a parameter-free comparison with the experimental data, 
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crossover to that of the gyration radius. It is required to con- 
sider the long time diffusivity (which has dynamic correla- 
tions) to explain the experimentally observed slow crossover 
of dynamic properties. Analytical calculations would be help- 
ful to obtain further insight into the results obtained here by 
simulations. 
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FIG. 3: Crossover of hydrodynamic radius: Comparison of theoret- 
ical swelling of static and dynamic measures (line with points) with 
experimental data (points) collated in Ref. 11511 . 

as shown in Figure for the measured values of qh against 
that of oig, leading us to the same conclusion. 
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FIG. 4: Parameter free comparison with measured values of qh 
against a g . Experiments are several data collated in Ref. fiUl . 

Based on the limited set of points shown in Figure |2j ef- 
fective exponents for the static and dynamic properties can 
be estimated at z — 10, which is typically the largest exper- 
imental value of z. Taking z oc \[M, we find R g ~ Af 59 , 
Ri ~ A/ - 58 , and R H ~ M 57 . These results should be taken 
only as indicative values of the exponent near z = 10, and are 
not conclusive critical exponents. 

To summarize, Brownian Dynamics simulations of the con- 
tinuum Edwards model were used to study the crossover be- 
havior of the hydrodynamic radius of a polymer chain. We 
found that the commonly used static measure of diffusivity 
(obtained from the inverse radius R{) has nearly identical 
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